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Abstract 

By means of the operator extension theory, we construct an explicitly solvable model 
of a simple-cubic three-dimensional regimented array of quantum dots in the presence 
of a uniform magnetic field. The spectral properties of the model are studied. It is 
proved that for each magnetic flux the band is the image of the spectrum of the tight- 
binding operator under an analytical transformation. In the case of rational magnetic 
flux the spectrum is described analytically. The flux-energy and angle-energy diagrams 
are obtained numerically. 

1 Introduction 

The problem of commensurability in periodic systems is a source of fascination in various 
branches of physics and is related to the appearance of the fractal spectrum for such systems. 
M. Ya. Azbel', D. R. Hofstadter, and G. H. Wannier PQ [H| were the first scientists who studied 
the simplest model of the two-dimensional (2D) electron system in a uniform magnetic field 
and discovered a self-similar energy spectrum, called "Hofstadter's butterfly". This picture 
shows the dependence of the energy on the magnetic field for two-dimensional periodic system 
in the tight-binding approximation and confirms the conjecture that the Hamiltonian of this 
model (which is reduced to the so-called Harper operator) has the Cantor spectrum at irrational 
value of the magnetic flux (the famous "Ten Martini Problem"). Very recently this problem 
has been solved affirmatively [3], [3]. 

It was assumed the long time that in the three-dimensional (3D) case the gaps in the 
energy spectrum, in general, are closed and the fractal structure of the spectrum is impossible. 
Recently M. Koshino, H. Aoki and co-authors have considered the dependence of the energy 
spectrum of a 3D periodic quantum system on the orientation of the magnetic field and on the 
strength of the field [S] [Bj- Using semiclassical methods they have reduced the Schrodinger 
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equation with the 3D tight-binding Hamiltonian to the ID Harper equation and shown that 
in the obtained angle-energy and flux-energy diagrams a fractal structure appears. Note the 
tight-binding method is based on a series of approximations of the initial Schrodinger equation 
that can change considerably the structure of the spectrum, e.g., the tight-binding Hamiltonian 
for a single-atomic lattice has only one magnetic band in the spectrum. 

Another explicitly solvable model of 3D periodic quantum systems was presented in [0] 
[TT] . We have considered spectral properties of the 3D Landau operator perturbed by periodic 
point potential. It was shown numerically in jTUj, JI] that the angle-energy and flux-energy 
diagrams for different bands of the spectrum have a fractal structure even in case of cubic 
lattice. Note that the dispersion equation in this model has a transcendental form and the 
analytical study of the spectrum becomes very difficult. 

In this paper we consider the quantum-mechanical model for the 3D periodic array of 
quantum dots obtained by means of operator extension theory CH]- It may be noted that 
this method was successfully applied for construction and spectral investigation of quantum 
periodic systems in two dimensions in our previous works JH], m particular, the model 
allows, in contrary to the tight-binding approximation method, to describe all bands of the 
spectrum. On the other hand, the dispersion equation of our model has more simple form 
than the corresponding one for the 3D Landau operator perturbed by periodic point potential, 
and in the case of rational magnetic flux the spectrum of 3D periodic array of quantum dots 
is described analytically. 

In the framework of our model, each Fock-Darwin level of a single quantum dot broads 
into magnetic band which in its turn splits into subbands. This is stipulated by the spectrum 
structure of irreducible representations of the magnetic translation group (for the first time it 
was shown by J. Zak PH])- If the numbers of the flux quanta of the magnetic field through 
each face of elementary cell of the Bravais lattice A of the system are rational (in this case the 
magnetic field is said to be rational with respect to the lattice A), then there is a finite number 
of subbands for a fixed magnetic band. If the fluxes through adjacent faces of elementary cell 
of the Bravais lattice of the system are equal to an irreducible fractions N±/Mx, N2/M2 and 
N 3 /M 3 , then the numbers of subbands are defined by the numbers Mi, M 2 and M 3 . Therefore 
the magnetic band arising from the fixed Fock-Darwin level can transform to a Cantor set 
when one of these fluxes tends to an irrational number. For the simple-cubic array, we prove 
that in both cases of rational and irrational flux the band is the image of the spectrum of the 
corresponding tight-binding operator under an analytic function. Hence, a fractal structure 
of angle-energy and flux-energy diagrams for 3D regimented array of quantum dots inherits 
from such diagrams for the tight-binding operator. It should be noted that similar description 
of the spectrum was obtained in the 2D case in , JH] and for quantum graphs in and 

Our model is directly relevant to the so-called 3D regimented arrays of semiconductor 
quantum dots which attract the increased attention lately [23] |2H]- The regimentation along 
all three directions results in the formation of an artificial crystal, where quantum dots play 
the role of atoms [23], [2H]- O. L. Lazarenkova and A. A. Balandin in [23] have analyzed the 
band structure of such systems using envelope function approximation. 
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2 The description of the model 



For construction of our model we start with the quantum-mechanical Hamiltonian H d of a 
single 3D quantum dot subjected to a uniform magnetic field B: 

H d = H + V(r), (1) 

where H is the Hamiltonian of a single electron in the magnetic field, 

Here e and m* are the charge and the effective mass of the electron respectively, and A(r) 
is the vector potential of the field B (i.e. B = rot A), V(r) is the confining potential of a 
three-dimensional quantum dot. In our model we choose V(r) as the potential of a spherical 
quantum well: 

V(r)=V(x,y,z) = ^f(x* + y* + z*). (3) 

Here the characteristic sizes d x , d y and d z of the quantum well are related to the frequency uj 
as follows: 

d x , d y , d z ~ >Jh/(m*uj). 

Now let us consider a periodic three-dimensional array of such quantum dots. We suppose 
that centers of the quantum dots coincide with nodes of simple cubic crystallic lattice A with 
the basis (ai,a 2 ,a 3 ) (Fig. 1). 

The state space of our model is the direct sum of the state spaces L 2 (IR 3 ) of single quantum 
dots: 

H = J2®H X , H X = L 2 (R 3 ) for each A G A. (4) 
AeA 

The following direct sum 

H° = J2 ®#a, H x = H d for each A G A, (5) 
AeA 

is the unperturbed Hamiltonian of the model. This operator is the Hamiltonian of a set of 
isolated quantum dots. To take into account the charge carrier tunnelling between dots we use 
the "restriction - extension" procedure of the operator extension theory |13j . We denote 
by T> the set of functions / from the domain T>(H d ) each of which vanishes in a neighborhood 
of zero. Let Sd be the restriction of the operator H d to the set V and let S is the following 
direct sum 

S = J2®S X , S x = S d for each A G A. (6) 
AeA 

We seek the "true" Hamiltonian H of the quantum dots array among non-trivial self-adjoint 
extensions of the operator S such as T>(H) n V(H°) = T>(S). 

The resolvent R of the operator H and the resolvent R° of the operator H° are related 
with the Krein resolvent formula (see, e.g. |15j): 

R(C) = R A (0 = R°(0 - g(C)[Q(0 + (C), (7) 

where operator valued holomorphic functions g(() and Q(C) are the so-called Krein T- and 
Q-functions of the operator H° respectively, and A is a self-adjoint operator in the deficiency 
space Z 2 (r) (the so-called "interaction" operator). The main part in the study of the spectrum 
is played by the Q— function Q(() and operator A. 
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2.1 The Q-function of the model 

In the considered case the Q-function Q(z) has the diagonal matrix in the space ^ 2 (A): 

Q(X,X';z) = 5 xy q(z), A, A' G A, (8) 

where q(z) is the Q-function of the operator Hd- To obtain the function q(z) we find prelimi- 
narily a general expression for the Q-function Q 3 (C', q\oo, B) of the operator at an arbitrary 
point q G M 3 . 

We denote by (u) the Hamiltonian of spherical quantum dot without magnetic field: 

H^\u;) = -^A + ^(4 + 4 + 4). (9) 

Let G 3 (q, q'; B) be the Green function of the operator H d = H^>(u, B), G 3 (q, q'; (\uj') be 
the Green function of the operator H^(u'), and Q 3 ((;q\u>') be the Q-function of the oper- 
ator HW(u'). We denote by S 3 (q, q'\u') and S 3 (q, q'\u, B) the singularities of the functions 
G^q, q'\ Cl^') an d G 3 (q, q'; (\co, B) at a point q respectively. It is known that for each to and 
lu' the singularities S 3 (q, q'\oj') and S 3 (q, q'\u, B) coincide and have the following form (see, 

e.g. my 

m* I 

S 3 (q, q'\u') = S 3 (q, q'\oj, B) = —- -. (10) 

2irn^ \q — q \ 

It implies for each uj' the following form of the function Q 3 ((; q\u, B): 

Q 3 (C; <?k B) = lim [G 3 (q, q '; - S 3 (q, q'\uj')} = 

q-^q 1 

= lim [G 3 (q, q>, C\u, B) - G 3 {q, q', C|a/)] + Q 3 ((; q\u'). (11) 

q-^q' 

We denote by H^ 2 '(u,B) the Hamiltonian of two-dimensional symmetric quantum dot sub- 
jected with magnetic field B (the field B is perpendicular to the quantum dot plane): 

H {2 \^ B) = ± (-ihV - e -A{r)) 2 + ^f(xj + x\). (12) 
and by H^ 2 \uj) the Hamiltonian of the same dot without magnetic field: 

H^\u J ) = -^A + ^f( X j + 4y (13) 

We use the standard notations uj c = \eB\/(cm*) for the cyclotron frequency and = ^/uj 2 + Auo 2 
for the hybrid frequency. It is easy to show that the eigenfunctions of the operators H^(u, B) 
and H^ 2 \VL/2) coincide and have the following form in the polar coordinates (p, 0): 



i&M) = ^d^i^TMy^ exp {~Q L!r ' (4) exp( ""^ 

m G Z,n = 0,1,2, (14) 

where Iq = ■s/KJ (m*f2), L^' is the Laguerre-Sonin polynomial of the power n. Corresponding 

eigenvalues EmniijJ, B) and E^(n/2) of the operators H^(cu,B) and H( 2 \tt/2) have the 
following form: 

B) = ^m+ ^(2n + \m\ + 1), (15) 



E%Km = E®{u,, B)-^m = f(2n + \m\ + 1). (16) 

Further we shall denote coordinates (p, 0, z) of a point q by (q±, q»), where q± = (p, (ft), q\\ = z. 
We have the following representation of the function Qs((; q\u, B): 



oo oo 
m=— oo n=0 



-Gfi(gn,ffii;C - Jsffi (n/2)|n/2)] + g 3 (C; g|n/2, w), (ir) 

where G\(x, x'\ (\<jS) is the Green function of the one-dimensional harmonic oscillator H^(u): 

„ m , . ft 2 d 2 m*cu 2 9 . . 

H^(u) = — + x 2 . 18 

v ; 2m* dx 2 2 v ; 

It is known that the function Gi(x, x'\ C\uS) has the following form (see, e.g. |27j): 
^ , i j-\ \ 1 x-, ( 1 ( \ tt I C \/2max(i, x' , , 

1 < ^> = 7a^ r U - s;) 17 (nb — ' x 



Therefore we obtain the following expression for the function Qs(C; q\oo, B): 

i / 2 \ 00 / 2 \ l m l 00 



(n + \m\)\ 



X 



Vis J 



i | 2(4S t (n/2)-Q 



xu „ (3*^,-*)) + Q3(C;,m.», (20, 

Now we turn to the desired Q-function q(z) = Q3{(; 0\u, B). It is known that for each u the 
function Qz{C,') 0|a;) has the following form (see, e.g. [26): 

m * — M 



Then it is easy to show that the g(z) has the following form: 

?(2) -^M^ r(i + £( B + i)-i)-?ife r(f + B - A )f + * (2) - (21) 
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2.2 The "interaction" operator of the model 



The matrix (A(X, A'))x,A'eA of the "interaction" operator A in our model is described by the 
following conditions (cf. [To]): 

1) A is the bounded self-adjoint operator in Z 2 (A), i.e. A(X, A') = A(X', A) for each A, A' G A; 

2) A is the nearest-neighbour-interaction operator, i.e. A(X, A') = 0, if |A — A'| > inf{|A — A'| : 
A, A' e A, A ^ A'}; 

3) A is invariant with respect to a natural unitary representation of the magnetic translation 
group in the space ^ 2 (A), i.e. 

A(X, A + n) = exp[7ri£ • (/x x X)}A(0, /x), A, /x G A, (22) 

where £ = B/§ is the vector of magnetic flux density (here $o — 27r/ic/e is the magnetic 
field quantum). 

Therefore in the case of simple cubic crystallic lattice A the matrix of the operator A is 
fully determined only by the elements A(X, 0), A G A: 



A(X,0) 



an, A = ±oi, 

a a , A = ±a 2 , (23) 

a 3 , A = ±a 3 , 

0, in the other cases. 



Here a±,a2 and a 3 are some real numbers (the so-called "coupling constants"). 

The construction of the Hamiltonian H = Ha of 3D periodic quantum dots array in a 
uniform magnetic field is complete. 



3 The structure of the spectrum of the Hamiltonian Ha 

In the framework of our model it is possible to show that the spectrum of the model Hamil- 
tonian Ha consists of the pure point and the band parts. The description of each band of 
the spectrum is mathematically reduced (without any additional simplifications or approxi- 
mations) to the investigation of the only energy band in the tight-binding model. Moreover 
in the case of so-called rational magnetic flux the structure of the bands may be described 
analytically. 

To investigate the spectrum of the operator Ha we start with the spectrum of the Hamilto- 
nian Hd of single quantum dot. The spectrum a (Ha) is discrete and consists of the eigenvalues 
E mn nx (the Fock-Darwin levels), 

U c ^ & 

E m n p n z = —m + — {2n p + \m\ + 1) + -{2n z + 1), 

mGZ, n p ,n z eNU{0}. (24) 

Thus, the spectrum of the Hamiltonian H° (the Hamiltonian of a set of isolated quantum dots) 
is pure point one and consists of the same points E mnpnz , and each point E mn nz is infinitely 
degenerated in the spectrum of H°. The Krein formula implies the following proposition. 

Proposition 1. Let a real number E be not a point of the spectrum of H°, then the following 
assertions are equivalent: 
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1) E e a(H A ); 

2) the operator q(E)I + A is not continuously invertible in the space l 2 (A). ■ 

Therefore, to describe the spectrum of the Hamiltonian Ha we can consider the spectrum 
of the operator A and the behavior of the function q(x) on the real axis. 
We denote by Htb the standard tight-binding Hamiltonian: 

Htb = -J2 exp(.i0\\>)CZC x „ (25) 

A, A' 

where the summation is over nearest-neighbor nodes of A, C x and C x are the standard creation 

e A ' 

and annihilation operators at the node A, 9 XX = — — f A-dl is the Peierls phase factor arising 

hc x 

from the magnetic flux. It is easy to show that the interaction operator A coincides with the 
operator H TB with t xx > = — A(X, A') for each A, A' G A. 

As for the function q, it is known that for real x the function q(x) monotone increases 
from — oo to +oo on each of the intervals (— oo,£o), (£q,£i), (^1,^2), • • • , where the numbers 
£0, £1, ■ ■ ■ are the eigenvalues E mripnz putting in ascending order (Fig. 2): 

£000 = £0 < £\ < ■ ■ ■ < £n < ■ ■ ■ ■ (26) 

Hence, there exists the multivalued real-analytical inverse function x{ x ) ( x £ ^) having 
continuos branches Xn(x) (n — 0,1, . . .) with the values in the intervals (—00, £0), (£0, £i),- ■ ■ (Fig. 3). 
Let us denote by r n the function r n (x) = —Xn( x )- 

The following theorem describes the structure of the spectrum of the Hamiltonian H A . 

Theorem 1. The spectrum a (Ha) of the Hamiltonian Ha is the union of two parts, ct(Ha) = 

£1 |J S 2 - The part £1 consists of infinitely degenerated eigenvalues E mnpTlz , where m^O. The 

00 

band part £ 2 of the spectrum of the operator Ha has the following form: £ 2 = U Y n , where 

n=0 

the set Y n is the image of the spectrum a(H TB ) of the tight-binding operator H TB under the 
real- analytical function r„ ■ 

Therefore the every band Y n is a result of deformation of a unique energy band in the 
tight-binding model under the corresponding function r n . This fact explains, in particular, 
various form of the spectrum in different energy bands. 

To continue the study of the spectrum of the operator Ha we consider the case of a rational 
field B. In this case the part £ 2 of the spectrum a(H A ) may be described analytically. 

Fix a basis (a±, a 2 , a 3 ) in A and denote 

B(a 2 x a 3 ) B(a 3 x ai) B(a ± x a 2 ) 

by the supposition, the numbers r]j are rational. Therefore, a basis (a[, a' 2 , a' 3 ) can be chosen 
in such a way that 771 = rj 2 = 0, 

B(a\ x a') N 

where N and M are coprime positive integers. 

Then by means of harmonic analysis on the magnetic translation group we reduce the 
spectral problem for the operator Ha to the eigenvalue problem of linear algebra. Namely, to 
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define the band part of the spectrum of Ha we have to solve the following dispersion equations: 



det[q(E) + A(p)] = 0, (27) 

where p = (pi,P2,P3) £ T 3 = [0,M _1 ) x [0,1) x [0,1) is the quasimomentum, A(p) is the 
M x M matrix with the following elements: 

A(m 1 ,m 2 ;p)= ^ A(X 1 a[ + (X 2 M + m x - m 2 )a' 2 + A 3 a' 3 , 0)x 

Ai,A 2 ,A 3 eZ 



-2rri (x x p! + \ 2 p 2 + X 3 p 3 + |Ai(MA 2 + m x + m 2 )) . (28 



exp 

It is obvious that for each p the series (J28J) is the finite sum. 

The following theorem describes the spectrum of the operator Ha in the rational case. 
Theorem 2. Let a basis (a^, a' 2 , a' 3 ) of the lattice A be choosen in such a way that 

B(a\ x at,) N , „. , 

V = V ^ 2) = m>^ B ( a 2 x a*) = B(a 3 x a[) = 0, 

where N/M is the irreducible rational fraction. Then the spectrum of the operator Ha is 
essential and is divided into two parts Si and E 2 . The first one, Ei, consists of infinitely 
degenerated Fock-Darwin levels E mn nz , where m/0. The second part E 2 is the band spectrum 
of Ha- This spectrum consists of bands I n , n > 0, lying on the intervals (e n ,e n+ i) (here 
e_i = — oo): 



E 2 — In-, In C (e n ,e n+ i). 



n=Q 



The each band I n , n > ; splits into M "magnetic" subbands I ni i, I = 1,2, ... , M: 

M 



n,l ; 

f=l 



where I n> i = {E n j(p) : p G T 3 } ; and /or /ixed p i/ie collection {E n>1 (p), . . . ,E n>M (p)} is the 
set of all solutions of the dispersion equation 

det[q(E)+A(p)] = 0, 

lying on the interval (e n ,e n+ i). 

The subbands I n j of the band I n can be overlapped. Some of subbands can be degenerated 
to a point ■ 

Therefore there are three kinds of points in the spectrum of the operator Ha- The following 
physical interpretation of the described structure of the energy spectrum of our model may 
be proposed. The points of the first kind (the levels of the single quantum dot) correspond to 
the classical orbits lying wholly in the single dot. The points of the second kind (the extended 
states) correspond to the classical propagation trajectories in the periodic system. Finally, 
the points of the third kind (the bound states satisfying the dispersion equation) correspond 
to the classical orbits lying in compact domain (cf. [Hj, |TT)]). 

As a result of numerical analysis of the dispersion equation (|27jl we obtain angle-energy 
and flux-energy diagrams represented on Fig. 4 and Fig. 5 respectively (we use a system of 
units such that c = h = e = 1 and m* = 1/2). According to Theorem 2, the bands of the 
spectrum look like deformed Hofstadter butterflies (cf. [T7|). 
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Figure Captions 

Fig. 1. Three-dimensional quantum dots array and basic vectors of array lattice. 
Fig. 2. The behaviour of the function q(x). 
Fig. 3. The behaviour of the functions r n (x). 

Fig. 4. The angle-energy diagram for the simple-cubic array of spherical quantum dots. 
Fig. 5. The flux-energy diagram for the simple-cubic array of spherical quantum dots. 
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